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Abstract 

We report on a high statistics simulation of SU (2) pure gauge field theory at 
finite temperature, using Symanzik action. We determine the critical cou- 
pling for the deconfinement phase transition on lattices up to 8 x 24, using 
Finite Size Scaling techniques. We find that the pattern of asymptotic scal- 
ing violation is essentially the same as the one observed with conventional, 
not improved action. On the other hand, the use of effective couplings de- 
fined in terms of plaquette expectation values shows a precocious scaling, 
with respect to an analogous analysis of data obtained by the use of Wilson 
action, which we interpret as an effect of improvement. 
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I. INTRODUCTION 



In the last few years the increased computer power available for lattice theorists has 
allowed a number of studies on relatively large lattices and it has been possible to test 
accurately scaling and asymptotic scaling. 

We have in mind in particular a series of works on SU (2) lattice gauge theory at finite 
temperature ( and references therein) which have been able to test for asymptotic 
scaling on lattices up to Nr = 16, showing still considerable scaling violations. 

On the other hand, in the same simulations one observes a rather good scaling of ratios 
of dimensionful quantities, thus supporting the idea that asymptotic scaling is violated 
by universal terms, which cancel in the ratio, and that some f3 function, other than the 
perturbative one, exists. 

Since the Renormalization Group, parameterized in terms of the "bare" coupling 
on the lattice, shows large deviations from the asymptotic behavior, we feel that it is 
very important to test different renormalization schemes, as stressed in a recent work 
by Lepage et al. 0. For instance, it has become more and more apparent that effective 
schemes on one hand improve the asymptotic scaling [^,0, and on the other hand reconcile 
the perturbative expansion of ultraviolet dominated quantities with the lattice numerical 
results 0. 

In this context, we have decided to test the effect of the use of a different lattice 
formulation, using an improved Symanzik action in the simulation of 4 dimensional lattice 
gauge theory. We did not expect a dramatic improvement of asymptotic scaling, even if 
in a pioneering work some indications of such an improvement were found, but we have 
considered that by comparison with similar works using Wilson action something about 
the origin of scaling violations could be understood. 

The study of this particular model, SU (2) gauge theory at finite temperature, is 
justified by its similarity with QCD, as well as by its greater simplicity. Moreover, the 
particular observable chosen, the temperature of deconfining phase transition, is well 
defined and can be exactly determined by use of Finite Size Scaling (FSS) techniques, 
thus allowing us to pursue our main objective, the study of scaling properties. 

We can anticipate our results: the pattern of asymptotic scaling violation observed 
with Symanzik action is essentially the same as the one observed with the Wilson for- 
mulation, the only apparent effect is a precocious transition from the strong to the weak 
coupling regime. This result seems to indicate that a universal lattice (3 function exists. 

Analogously the effective coupling approach, formulated in terms of plaquettes ex- 
pectation values, works well in improving the asymptotic scaling figure and gives results 
consistent for the two lattice formulations. 

In section H we give a short review of the Symanzik approach to the lattice formulation 
of field theories and we recall physical characteristics of the model under study. 

In section |T|, we review the FSS techniques used in this work, in order to give a 
presentation as self contained as possible. 

In section |V| we give details about our simulations. 

In section we present the results of our measurements, exploiting the FSS techniques 
to present various tests of consistency of the results. 

We discuss the results in section together with an analysis of the different renor- 
malization schemes. 



Appendix A is devoted to technical details on the subtraction of biases and on error 
estimation. 

Analogously in Appendix B we discuss our implementation of the Density of States 
Method (DSM). 

II. THE SYMANZIK IMPROVEMENT AND SU (2) LATTICE GAUGE 

THEORY 

1. Generalities 

The formulation of a field theory on a lattice is a necessary step if we want to study 
it in a non perturbative way with a computer simulation. The lattice can be seen as an 
ultraviolet regulator as any other, introducing in a natural way a short distance cut-off. 

We can extract from a computer simulation relations between observables (for ex- 
ample, mass ratios) only for this regularized theory, and obviously we cannot recover 
completely the continuum limit, so our predictions will be affected by systematic errors 
(the so-called lattice artifacts). 

The cut-off dependence of a theory can be reduced by a clever choice of the regu- 
larization scheme: this is a well known technical point which was studied extensively by 
K. Symanzik He realized that this fact can be used to minimize the consequences 

of a non-zero lattice spacing a (for an introductory discussion see the pedagogical work 
of Parisi g). 

The key observation is that a field theory can be transcribed on a lattice with consid- 
erable arbitrariness: the lattice action must only reduce to the continuum one as a 0. 
So we can redefine the action by adding an arbitrary combination of irrelevant operators, 
which vanish in this limit, and this is equivalent to change the regularization scheme. 
Symanzik has also shown (see for instance [0) that every lattice regulated theory is per- 
turbatively equivalent to a local effective lagrangian, in a given renormalization scheme, 
of the type 

.Ciat. ^C,. = T. +a^Y. ^Of^ + E <^Of^ + ..., (1) 

ill 

where ol"'^ are local operators of dimension d + n, if is the dimension of space. It follows 
that it is possible to build a lattice action in such a way that the corresponding effective 
lagrangian is 
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using an appropriate linear combination of the operators o'f^'^^^ transcribed in an arbi- 
trary way on the lattice. The determination of the coefficients for this combination is 
a non trivial task for an interacting theory: in principle it could be done with an high 
precision numerical simulation. 

We can also compute them perturbatively, if the bare coupling constant is sufficiently 
small (i.e. near the continuum limit for an asymptotically free theory), using a matching 
procedure between vertex functions 0, or in the case of gauge theories (for instance) 



matching the expectation value of gauge invariant quantities, hke the Wilson loops 



or the inter-quark potential [jT2| . A exposition of the different improvement strategies can 
be found in the work of Liischer et al. ||13|| . 

If this program works, we may hope to obtain more accurate predictions for physical 
quantities with the not extremely big lattices we can use with the current computing 
resources. The more direct way to verify that lattice artifacts are reduced is to measure 
an adimensional ratio of physical quantities while making a smaller and smaller. For 
example using the action in Eq. we expect to find for a mass ratio = correlation 
length) 
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(3) 



A more refined test consists in verifying how well the quantities A4i we can measure on 
lattice obey the Renormalization Group equation, which is merely a statement of cut-off 
independence of physical observables (scaling). The lattice artifacts modify this equation 
adding non-universal (i.e. A^-dependent) scaling-violating terms, so we have 
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Here (3m is a non analytic contribution that cannot be calculated perturbatively (for 
instance, in the bidimensional a model it can be evaluated in the large N limit, where 



one finds (3m = O 



and A 
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If we use the action 



with 



p = 2, we expect first of all Am = O (a/C,)^ In {a/^) , and then a reduced Pm- In the real 
case we can calculate the improved action at best to the first perturbative order, and we 
expect contributions O (7^ (a/^)^ In (a/^) to A^, which give a little contribution for an 
asymptotically free theory if we are sufficiently near the continuum. 



2. Scaling vs. asymptotic scaling. 

In the general case we know only the first few perturbative terms of the function 
^{g), in particular the first two (scheme independent) terms. These are the relevant 
terms in the continuum limit if Asymptotic Freedom holds, while in an intermediate 
coupling range higher order contributions can be important. If an observable Ai follows 
the renormalization group with the perturbative approximation to (3{g) we say we are in 
an asymptotic scaling region: if this is not true, it does not imply that the data cannot be 
trusted, but only that we had to use a better approximation to the exact beta function. 
We can try to extract this "improved" (3{g) from the simulation itself and then verify 
the self consistency of this procedure. If we can find a function Pesig) so that the RG 
evolution holds for different quantities we say we are in the scaling regime. 

It is important to note that 

• The Symanzik program can improve the scaling, but not necessarily the asymp- 
totic scaling: indeed, the onset of asymptotic scaling is expected in a region where 
irrelevant operators give negligible contribution. 



• The scaling test is more conclusive than the stabilization of mass ratios, as we can 
imagine that there are some Ai independent contributions to Am which cancel in 
Eq. (I). 

As it is well known asymptotic scaling in the SU (2) gauge theory with the Wilson 
action has not yet been found (see for instance So we must extrapolate the beta 
function in some physically sensible way in the intermediate coupling region. 



3. The model under study 

The SU (2) gauge theory at finite temperature has been discussed at length for in- 
stance in |T^: we here just recall that the model shares with QCD the presence of a 
deconfining phase transition, characterized by the breaking of the center Z (2) symme- 
try |jl5| and by the appearance of a non zero expectation value for the thermal Wilson 
loop, or Polyakov line, P. It is by now well established that the transition is of the sec- 
ond order and all numerical simulations give results in good agreement with the early 
ansatz [|l^ that the model belongs to the same universality class as the 3D Ising model. 

In our simulations we have used a "tree improved" action, i.e. we have only corrected 
in part lattice artifacts in the classical theory. 

This can be done by observing that the usual Wilson action is equivalent to the 
effective lagrangian 



= -\f^uF^, + ^a^d^F^^d^F^, + O(a^). (5) 



We can compensate for the O(a^) term by adding a suitable irrelevant operator. A possible 
and widely used choice is 

Sj = Sw + = I3Y^U,^,+(3J2 Quixi - ^f^ixs) = /9E {l^ixi - ^^^1x2) , (6) 

where Unxm is the n x m plaquette. 

This choice has already been experimented in other works, both in the study of the 
finite temperature theory |^ and of the quark-antiquark potential , using the ichosae- 
dral approximation to the SU (2) group. To our knowledge the present work is the first 
where a 4-dimensional gauge theory is simulated with the Symanzik action and the full 
group. 



III. FINITE SIZE SCALING ANALYSIS 



The Finite Size Scaling (FSS) technique is by now a widely used tool in the investi- 
gation of pseudo-critical properties of statistical systems in finite volume. Its application 
to the analysis of second order phase transitions permitted a detailed numerical test of 
theoretical predictions on critical exponents. The effectiveness in the study of SU{2) de- 
confining phase transition has been already demonstrated in a number of papers, where 
use of the Wilson's action was made |I|,p|JT8[]. 



In this section we collect some formulas needed by our investigation, in order to 
give a presentation as self contained as possible of our results. A more detailed general 
introduction can be found for instance in |T^, while we refer to the original works 



for the improvements of the method which we have applied to our analysis. 

Numerical simulations of statistical systems are limited to finite volumes, character- 
ized by some length scale L. In the vicinity of an (infinite volume) phase transition 
the system exhibits a "pseudo-critical" behavior, for instance at a second order phase 
transition the susceptibility x shows a peak broadened by the finite volume, while the 
correlation length ^, defined in terms of appropriate pair correlation functions, reaches 
the dimensions of the system. For our purposes, the shift in the critical temperature is 
the most important effect, which can be easily estimated if one assumes some universality 
class for the model under consideration and that the correlation length exhibits a power 
behavior in the vicinity of the pseudo-transition. 

i^\tr (7) 

where t = is the reduced temperature ( in units of the infinite volume critical 

temperature ). If the finite volume "critical" temperature, Tc{L), is reached when the 
correlation length is of the order of the physical size of the system, it follows that the 
expected shift is 

|T,(L)-T,| cxT.L^^ . (8) 

Note that in fitting measurements the critical temperature is a very sensitive parameter, 
hence the error in its determination should be reduced as much as possible if one is 
interested in determining universal parameters with a good accuracy. 

To extract infinite volume limits it is necessary an "ansatz" on the form of the ob- 
servables. 

The main idea is that near a second-order phase transition, as the correlation length 
^ approaches the size L of the system, the observables O depend on this two scales only 
through the ratio r = L/^, apart from a volume dependent prefactor 

Ol = L'^Qo (r) (9) 

or, in terms of the reduced temperature, 

Ol = L^Qo {Lt") . (10) 

If for large L this observable exhibits a critical behavior, with exponent —p, the 
function Qo should behave like some power of its argument, and the requirement of 
cancellation of volume dependence requires that u = p/u. Hence it results the following 
parameterization 

Ol = LpZ-^Ql {Lf) , (11) 
or alternatively, and more commonly 



Ol = Lp/^Ql {tL"-) 



(12) 



The important point is that renormahzation group arguments show that all the critical 
behavior can be derived from the singular part of the free energy density, which is then 
assumed to have the conventional form 



/, = L-'Q (grL'/'^, g^L^) (13) 



where qt, Qh are connected to the reduced temperature and the external magnetic field 
by linear relations, plus corrections 



gt = Ctt + [th, 

gh = Chh + 0(th,h^) . (14) 



The additional dependence on irrelevant scaling fields gi, with negative exponents yi, 
cannot be neglected in many practical simulations, and determines corrections to the 
scaling behavior. Note also that we name t] the critical exponent of the "magnetization", 
instead of the conventional (3, to avoid confusion with the coupling in the latticized theory. 

Let us come to our system, and recall that we simulate on a lattice of volume L^, with 
L = N„a, and temporal extent Nt, determining a temperature T = 1/ (Nta). 

The relation between the spacing a and the bare coupling /3 = is unknown, and 
one of the motivations of our simulation is the study of the validity of the universal 
two-loop asymptotic formula 

''^-(24)'"'^'-p(-44)- 

As first noted in |T|, it is convenient to rewrite the free energy density in terms of the 
dimensionless combination LT, that is, of the ratio y = Nfj/Nr. 

^{t,h■N^■Nr) = y-''Qf{gtyKghy'^,9^y''') • (16) 
Define P the spatial average of the Polyakov line 

Nt 

L{x) =Y[U (f , t) 
t=i 

P = K%^^[L{x)] (17) 

and introduce a source h for this quantity in the partition function, through the operator 
hZ (a, Nr) N^P, where the renormahzation factor Z cancels the divergent self-energy 
contributions to the Polyakov loop. 

Derivatives with respect to the source h allow to define the physical order parameter^ 
and the susceptibilities 



In a finite volume system, as tunnelling between inequivalent "vacua" tends to restore the 
Z (2) symmetry, the experimental order parameter is defined as \P\ 



(P) 
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where use has been made of the first hyperscahng relation 

7/1^ + 2rj/v = d 



(19) 



implied by Eq. (|T8D. Q 

A useful quantity in the determination of the critical temperature is then the Binder 
cumulant (74, defined as the fourth derivative of the free energy with respect to "magnetic 
field" , normalized to the susceptibility 



94 



d'fs 



h=0 



(20) 



In fact, its expression is a directly scaling function, where explicit dependence on N^r is 
canceled by the use of the hyperscaling relation: moreover, the multiplicative renormal- 
ization Z of the Polyakov loop cancels, and therefore the expression for (74 can be safely 
given in terms of expectation values of moments of Polyakov loop 
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(21) 



Its directly scaling form can be expressed as 



2 



In fact, assuming that fluctuations contribute to the singular part of the free energy density 



as the ratio of the unit volume to the volume ^ {t)'^ 
which is incorporated in Eq. ([TBI), if for small values of x 
Now derivatives in Eq. (|l8|) can be written as 

(p) = - (ivjivo"'^^"^^^/'^ Qi {gt {NJNrf-' 



and so, combining the dependence of Qi^2 on gt (that is, the same as Q/ ), with the requirement 
that in the infinite volume — > 00 limit the critical behavior is specified by the exponents r/, 7 

it results Eq. (0) . 



9i 



it; Nr) = Qg {gt {t, N^) y'/", gw^') (22) 



where we have conserved the first irrelevant scaling field. The localization of the crit- 
ical point is then possible by finding the intersections of (74 curves, as functions of the 
temperature t, at various values of the spatial size N„. 

Taking into account the first irrelevant field one expands the expression given in 
Eq. ( p2D in the vicinity of the transition 

g^ {t; N^; Nr) ~ ^74,0 + g4,ity^ + ^74,2?/^^ • (23) 

Hence the intersection of two g4 curves at values of the spatial size N„, N'^ corresponds to 
a shifted temperature 



94,1 



where b = N'^/N„. 

As ?/i < 0, this shift in the crossing goes to zero as N^j grows, while also at fixed A'^^- 
as b grows the shift goes down, governed by the exponent — 1/z/. 

Actually, in the simulation of SU {N) gauge theories, the connection between the bare 
coupling f3 and the lattice spacing, which allows to determine the temperature in units of 
the reference scale A, is known only in the asymptotic regime go ^ 0. As noted in this 
means that expressing the reduced temperature t in terms of /3 through an approximate 
formula like 



^Cabo 



2Cab, , 



(25) 



introduces an error O in the determination of t. This shift is of the order of 8% in 

the relevant coupling regime explored, but one should be aware that this estimate would 
be a posteriori justified by the observation of an asymptotic scaling behavior, which is 
absent. 

In our search for the determination of the critical temperature this shift is immaterial, 
as long as a relation of the form 

t oc (/3 - /3,,oo) (26) 

is valid. 

Therefore Eq. ( ^4|) translates in the form 

(3,{N„,K)=(3,,^{l-ce) (27) 

where 

and the constant c maintains a dependence on A^ and /?c,oo which is irrelevant in the 
£ — i> extrapolation. 



IV. DETAILS OF THE SIMULATION 



Most of the work has been done on the APE supercomputer ||20|. The model operating 



at Pisa is the so called "tube" machine, a 128 processor parallel computer with a peak 
performance of 6 GigaFlops. 

All the code has been written in the high-level language proper of this machine 
( "apese" ), and the program runs at about 35% of the peak speed: this result is quite good 
if one takes into account the complications both in addressing and in memory access due to 
the use of the Symanzik action. In other words, in our code the floating point performance 
is slowed by the large number of integer and addressing operations required to load the 
elements of group in the register- file: this means that an analogous implementation for the 
SU (3) gauge group should result in a better performance, as the ratio of floating point 
to integer operations would be increased both in computation and in I/O. A detailed 
account of the implementation will be given elsewhere. 

In the update of lattices with A^,- = 3, 4, ... 8 we have used an overrelaxed Heat Bath 
update. An exact Heat Bath algorithm has been implemented, in the modified Kennedy- 



Pendleton form which results in higher acceptance and is well suited to a parallel 
machine like APE; a number of complete overrelaxation sweeps ranging from 10 to 16 has 
been used to decorrelate between subsequent Heat Bath sweeps. 

As the Symanzik action couples next-to-nearest neighbors sites, on a lattice 2 x Ncr 
it results a self coupling of spatial links, and therefore the action is quadratic and pre- 
vents the use of an Heat Bath algorithm: hence on these smaller lattices we have used a 
Metropolis algorithm running on RISC workstations. 

In Tabb. |,fT|,. . . |V| we present the simulation parameters, together with the estimated 
autocorrelation time of the Polyakov line. The Density of States Method has been used 
to interpolate between simulated data points. See Appendix B for a discussion of the 
meaning of the last column, where the "reweighting range" used in the interpolation is 
listed. 

We have performed also a series of runs on large symmetric lattices to measure the 
expectation values of plaquettes, Uixi, t/ix2, to be used in determining effective couplings 
"a la Parisi". We list in Tab. |V11| the corresponding run parameters. 

V. RESULTS 

A. Determination of critical point by FSS techniques 

Runs on lattices with A^^ = 2, 3, ... 8 and up to N„ = 24 have been devoted to explore 
the transition region and measure accurately the Polyakov line. 

1. Binder method 



In Figg. |l|Jl,. . . 1^ we show the plots of the Binder cumulant, Eq. (21) for different 
values of Nr and N„. 



As the Binder cumulant is not a self-averaging quantity, the resulting bias has been 
computed, and found to be relevant for the larger lattices, where statistic is relatively 
poor: details on the evaluation and subtraction of the bias are reported in Appendix A. 

The solid lines, as well as the dashed lines which are an estimate of the error, are 
obtained with our implementation of the Density of States Method P5|P^ , described in 
Appendix B. 

The determination of the crossing points as well is based on DSM and we report in 



Tab. [VIII| the values of intersections /3c {N^, Nfji) for the different lattices. 

We make use of the extrapolation in Eq. (|27|) based on the introduction of a single 
irrelevant field, with exponent y = —1, to obtain for lattices 3 — 6 the infinite volume 
limit of the critical coupling: 

Nr = 3 Pcoo = 1.59624(13) a = 0.23(5) 

Nr = 4: = 1.699(1) a = 0.20(6) 

Nr = 5 f3c,oo = 1.76948(3) a = 0.85(10) ^ ' 

Nr = Q /3c,oo = 1.8287(11) a = 3(1) . 

For what concerns the lattice with A^,- = 2, we have not enough statistics to distinguish 
the different intersection points. As this point is in the strong coupling region, an high 
precision determination is not necessary. 

For the Nr = 7, 8 lattices we have only a single crossing, so we can estimate the 
critical coupling by assuming some value for the a coefficient. By choosing a ~ 3 — 5 for 
the Nr = 7 lattice we obtain, having b = 20/14 = 1.43, a value for e ~ 4.2 x 10~^, hence 
ae ^ 1 — 2 X 10~^. Analogously one obtains ae ~ 3 x 10~^ for the N^- = 8 lattice, so we 
estimate a systematic error of order 10~^, giving so 



A^^ = 7 /3e,oo = 1.8747 ± 0.002 ± 0.002 
A^^ = 8 /3c'oo = 1.920 ±0.004 ±0.003 . 



(30) 



Using the two loop asymptotic formula, Eq. (|^), this values are reexpressed in terms 
of the T/Ai ratio in Tab. IX . 



2. X method 

A different way to determine the critical point has been recently proposed by Engels et 
al. in [|18|. The idea is to exploit the specific form of the scaling law for the susceptibihty, 
which can be written as follows 

X. = Nl{P^) , (31) 

taking into account that in finite volume the expectation value of Polyakov loops is set to 
zero by spin flips between degenerate vacua. 

Expanding in the vicinity of the transition the expression in Eq. (0) one obtains 

Xv= \ ^] ^ Co ± Ci ± C2 — t— ± C3 — 



ln(x.)U^ln(co) + 2.n(^)+5j(^)"'. (32) 



Hence, assumed that the effect of the irrelevant is small [yi ~ —1), at the phase transition 
the true susceptibility Xv, plotted as a function of the spatial size on a doubly logarithmic 
scale, should exhibit a linear behavior. 

Let us apply this method to the Nr = 3 lattice. We use the density of states method 
to determine Xv at various values of /3, and fit the data to a doubly logarithmic law. The 
resulting value for the minimum of is at f3c = 1.5956(4) at 95% c.l., which is in good 
agreement with the determination. One obtains also the slope 7/z/ = 1.91(3), which 
exhibits a 3% difference from the known value in the 3D Ising model ^jv = 1.965(11). 

The same procedure has been applied to the other lattices: we find that, when appli- 
cable to our data, this method gives essentially the same result as the method of Binder 
cumulant, so we have used it as a cross check of the validity of the analysis. 



3. Universal Scaling Behavior 

We may test the universal scaling laws, expressed by Eqq. (122|) and (H), by plotting 
g4 {t, y), { y = N^/Nr ), as a function of the combination ty^^" , and the analogous expres- 
sion for the Polyakov loop, N^-y'^^'^ we choose u = 0.628, rf/i' = 0.516, 7/1^ = 1.965 
as given by the 3D Ising model. In Fig. ^ one may see that the values cluster round a 
straight line, thus supporting the theoretical conclusions. The agreement is not so good 
for the Polyakov line. Fig. ^, which is sensitive to the self energy contribution. 



B. Critical temperature in the bare scheme 



Let us give in Fig. |T0| a cumulative plot of the various determinations of the critical 
temperature. The data on Wilson action are taken from |l[] and normalized to the point 
at A^^ = 8, ( instead of using, for instance, the perturbative A ratio). 

A few comments are in order: 

• the weak coupling regime appears to be reached earlier in the simulation with the 
Symanzik action, 

• in the weak coupling regime the pattern of asymptotic scaling violation appears to 
be essentially the same, by using Wilson and Symanzik action. 

This result is in a clear disagreement with the earlier work by Curci et al. 0, where 
asymptotic scaling was found with the Symanzik action already at A^^ = 5. We have then 
analyzed our data with the same method, based on a fit to the critical behavior of the 
Polyakov line. For (3 > (3c, the form 

{\L\) = A iP - (33) 

is assumed, while for (3 < (3c the Polyakov line is set to zero: this introduces some arbi- 
trariness in the fitting procedure. In fact, using Eq. ( P^D to determine (3c one is forced to 
discard data points assumed to be on the "left" of the transition. Moreover, on the "right" 
side of the transition (smaller lattice spacing) there are stronger renormalization effects, 
and these introduce a systematic which could be reduced by subtracting perturbative 
tails. 



In Fig. ^ we show the data points, together with the results of a three parameters 
fit to the critical behavior: the study is done for lattices 3 x 10,4 x 12,5 x 16. 

Results of the fit are shown in Tab. together with the value of the best reduced 
found. Errors are estimated with the method of equisurfaces. 

The resulting critical temperatures are also reported in Fig. [1^. The high value of 
the reduced x^; especially on the 4 x 12 lattice, shows the bad quality of the fit, but for 
our purposes it is important to point out that the finite volume "critical" temperature 
determined by this method shows no sign of asymptotic scaling. On the contrary, it is 
present only a constant shift from the infinite volume limit, and we conclude that the 
results of were hampered by a insufficient statistic. 



1. Perturbative A ratio 



We can exploit the high precision data obtained with the Wilson action by Fingberg et 
al. |l[] to give in Tab. ^ the resulting determination of A ratios, by using the asymptotic 
formula 

These values are to be confronted with the perturbative results from Weisz and 
Wohlert ^ 

^ = 4.13089(1) . (35) 
Al 

showing a discrepancy of the order of 10%. 



VI. DISCUSSION OF RESULTS 

We can resume the outcome of our simulations as follows: Symanzik and Wilson 
actions give consistent results in the so called "bare" scheme, that is, the Renormalization 
Group evolution driven by scheme independent part of the perturbative f5 function does 
not seem to show scaling of the temperature. The pattern of scaling violation is similar, 
at least for iV^ > 4. 

On the other hand it is well known that adimensional ratios of physical quantities show 
a good scaling in the coupling range considered in this work, at least in the case of Wilson 
action, where the ratio Tc/a^^"^ has been showed to scale over the entire P = 2.30 — 2.74 
range |l|. 

In our opinion these two facts, that is, the good scaling and the similar pattern of 
asymptotic scaling violation in the two actions, are consistent: as anticipated, in the 
scaling region the improvement has no effect on the asymptotic scaling. 

Let us first elaborate on the statement that the two scaling figures are compatible 
with each other: to do this, we set up a perturbative framework and discuss the effect of 
the coupling redefinition implied by the use of Symanzik action. 



We split l3{g) in the universal = —h^g^ — big^ and in the scheme dependent (un- 
known) part (3s.d. = 0{g'^), where 

(47r)' 3 (47r)' 3 ^ ^ 

The beta function determines the dependence of lattice spacing a on the g value: 



<S) = ^ (M") exp I ( 1 + 9^ <i„/" I . (37) 



bi / 1 ^ ^ °° 

. , , ^'o exp ; — 

The coefficients dn in Eq. (|37D are generated by Z^^.^.; they play the role of asymptotic 
scaling violating terms. For instance, in terms of the first sub-asymptotic term in f3s.d., 
the 62 coefficient, one easily finds 

If we redefine the coupling constant 

9 = g{g) =9 + (yif + oi^t + ■ ■ ■ + + • • • (39) 

we obtain 

<9) = (bof) ^ exp -— ^ ](^+9"i: dnf' ■ (40) 



n=0 



Al ^ ^ V 26o^' 

This is the evolution dictated by a new beta function P*{g), 

PTg) = ma)) = m) + KAa) , (4i) 

and the non universal terms are changed, for instance one has 

6* = 62 - 2ai6i - 3a?6o + 26o«2 • (42) 
The prefactor exp (ai/60) redefines the scale parameter, consequently 

ai = -ho\n(^\ . (43) 



vAl/ 

On the other hand the asymptotic scaling violating contribution is changed 



and one may hope to reduce its effect by a clever (and physically motivated) choice of g. 

To study if the asymptotic scaling violation can be ascribed to the d2 terms, we have 
performed an exploration of values of 62 for both actions: in Fig. ^ we show how setting 
values of 62 for Wilson and Symanzik actions respectively to 0.0018, 0.0011, one is able to 
obtain a good scaling: the N^- = 2 term has been dropped for clarity. 

On the other hand, we are well aware of the arbitrariness of the procedure, as we have 
checked that this fit procedure is unstable if also the coefficient 63 is introduced. This is 
easily understood as values of 62 so large are hard to justify on a perturbative ground: 
too see this consider the running coupling equation in terms of the characteristic QCD 
coupling, as/ (47r) 

--f?^)^-f?^)^4330 f^V + ... : (45) 

for instance at /? ~ 1.9 one has as/ (47r) ^ 0.013, which means that the last term in 
Eq. (|45| ) is of the same order of magnitude as the preceding one. So we can merely regard 
this procedure as the definition of an effective 62, much in the spirit of the effective A used 

in m. 



2. Perturbative approach: q„ coupling 

Let us consider more closely the poor convergence of the speculative expansion given 
in Eq. (^5]): in the preceding subsection we have looked for an appropriate 62 coefficient 
without any physical motivation. On the other hand, Lepage and Mackenzie point out 
how a bad choice of the expansion parameter may give rise to a poorer convergence of the 
perturbation series, and they suggest strategies to design improved couplings, defined in 
terms of physical quantities. 

This is much in the spirit of the old suggestions of Parisi, to define an effective coupling 
by inverting the perturbative expression for some ultraviolet-dominated quantity, like the 
single plaquette. 

In [0] it is suggested to define an effective coupling a^ in terms of the quark-quark 
potential, by the relation 

F(g) = -^^i^%^. (46) 

This coupling can be perturbatively related to ao, the bare lattice coupling, by the 
formula 



a, 



if) = ao 1^1 + ao 1^(477) h In i^^^ + ^) ) ' ^^'^^ 



where the constant K can be extracted by Kovacs work [E^]: she reports 



Numerically one finds 



Jr 



-16.954, for = 2 
-19.695, for = 3 



(49) 



To translate the result for the Symanzik action it is sufficient to rescale A factors and 
impose equality in the physical result. The defining relation is 



ho In 



(47r) 



A?, 



(47r 



One obtains easily, for SU{2), 



11 



^/ = ^L + ^ln(^ 



(50) 



(51) 



It is then possible to use the result of Weisz and Wohlert [25| on the A ratio: 



A/ 



exp 



(47r)6o 
4.13089(1), for A^ = 2 
5.29210(1), for A^ = 3 



Ca (0.04329017(1)) - (0.04141417(1)) 



and it follows 



Jr 



-6.54859, for A^ 
-7.47609, for A^ 



2 
3 



(52) 



(53) 



Translating into the expression for K, one easily obtains 

Kl- 



2.69831, for A^ = 2 
4.7018, for A^ = 3 



1.04224, for A^ = 2 
1.78479, for A^ = 3 



(54) 



We can now use value of ao to set the value a„ (vr/a): one has obviously 

(vr/a) = ao (l + tto-K'LC/)) • 

For a given value Pc = Ca/ (27rac) of the measured critical coupling, the corresponding 
critical temperature in the ay scheme may be obtained by using the relation 

T _ 1 vr 
A„ vrA^^ aAy 

and the two loop scaling formula 



(55) 



(56) 



- = exp 

1 Vy 



— TTT In 

26^ 



^Airboay (q) ^ 



(57) 



^vr) bo ay (q) ^uo 

Let us report in Tabb. pCII| and [XIII| the results of the analysis both for the Wilson 
and the Symanzik case. 

It appears that the use of a perturbative expression for the effective coupling ay does 
not give a better asymptotic scaling. It should be stressed however that this was merely 
an exercise, as the Lepage scheme is intended to work in conjunction to a non perturbative 
determination of the coupling ay, based on the measurement of the heavy quark potential. 



3. Non-perturbative approach 



The bare coupling itself has no meaning in itself, so it must be eliminated in favor 
of some physical quantity, connected to it. In the perturbative approach this is a way to 
hide infinities and obtain a renormalized field theory. 

A simple way to implement this program on the lattice is the elimination of the 
bare coupling in favor of some UV dominated quantity, such as the expectation value of 
plaquettes. To this end, let us recall the results of the work of Weisz and Wohlert [^] on 
their perturbative expansion: by defining, for simplicity in the adjoint representation 

1 -| cxD g2n 



In 

the first coefficient is given by 



(Tr [/(L,T)) 



- i:^-»(L,r) , (58) 



Wi{L,T) = C,I{L,T) (59) 
and the following values are given 

^ ' ^ \ 0.366262 Symanzik 
rf. .^ ^ / 0.862251 Wilson 

^ ' ' \ 0.662624 Symanzik ' ^ ^ 

Coefficients for the Wilson case are known up to fourth W4, while for Symanzik we shall 
limit ourselves to lowest order. 

This allows to define two effective couplings in the Symanzik case, depending on the 
two expectation values measured. That is, knowing that 



In 



1 1 



— (Tr[/(L,r)) 



jCfI{L,T) (61) 



and given /3 = we define 



Pl(2) - <~-aCf- . (62) 

-L — L/ixl(lx2) 

Numerically, for the SU{2) case, this corresponds to 

0.549393(0.993936) 
A(2) = — - — Yi ■ 

-L — t;ixl(lx2) 

By using the results of measurements on large symmetric lattices we can obtain via 
interpolation the values of plaquettes needed to compute the effective coupling: we report 
the results in Tab. ( pClV| ). 

We plot in Fig. |13| the comparison of the different schemes, together with the Wilson 
data, taken from [0], in the effective scheme deduced from the single plaquette. As usual 
we normalize to the last point: we drop the point at A^^ = 2 to show in a larger scale the 
other points. It is worth noting that 



• the two effective schemes used for the Symanzik action agree with each other, 

• they both show a good scahng starting at A^i- = 4, 

• the same scheme used with Wilson action seems to show, starting at A^,- = 5 the 
same behavior. 

We think that these results can be interpreted as a confirmation of the value of effective 
schemes, and as an indication that Symanzik action gives an enlargement of the scaling 
window. 

VII. CONCLUSIONS 

In this work we have performed an high statistics simulation of the SU{2) pure gauge 
theory at finite temperature, using the tree-improved Symanzik action and exploring a 
range of lattices up to A^^ = 8; we have determined, with the help of Finite Size Scaling 
methods, the critical couplings for the deconfinement transition, with the purpose of 
studying the scaling properties of the critical temperature and the effect of improvement. 

We can resume as follows our findings: 

• asymptotic scaling violations are present, and are of the same size as in the case of 
Wilson action. The pattern of violation is similar and can be interpreted as driven 
by the same lattice P function, modulo scheme redefinitions. This means either that 
lattice artifacts are small, or they are the same for the two actions, in accordance 
with previous determinations of the scaling window. 

• The use of a non perturbative coupling derived from the plaquette expectation value 

gives a better scaling figure: this effective coupling was already used in the same 
context but with the Wilson action, and comparison of the two analyses shows that 
Symanzik action apparently gives a slightly precocious scaling. 

• In our study we are unable to directly test the scaling, so it could be desirable a zero 
temperature study of some other dimensionful quantity, like masses or the string 
tension. 

From a practical point of view, if we are willing to trust the precocious onset of 
scaling in the effective scheme, we can say that the gain in volume obtained by use of 
the improved action compensates the increased complexity in the update routine, thus 
slightly favouring the use of this technique in simulations. This should be even more true 
in the case oi SU (3) group. 
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APPENDIX A: STATISTICAL EVALUATION OF BINDER'S CUMULANT. 



In order to extract in a reliable way the quantities we are interested in from a set of 
N measures we have to correct for the distortions induced by the finite sampling. It is 
well known that often the "naive" estimator for a statistical quantity is not the best one. 
The much simpler and well known example is the estimator of variance cr^ for a set of N 
uncorrelated data Xj. One may try 

E,{a') = >) - £^(< X >)^ = (1 E - E E I (Al) 

but, owing to the fact that E{< x >) is a random variable with a finite variance the 
estimator's expectation value is 

T N 1 AT / 1 \ 

< M^') >= ^E < > . E (% <^'> +(1 - 5.,) <x = j^) a\ 

% — 1 i — I5J — 1 

(A2) 



which is "biased" by a finite sample effect 0{j^). To evaluate in the best way cr^ we must 
correct for this bias: this can be done using a new redefined estimator 

E{a^) = N/{N-l)E,{a^) (A3) 

which gives the correct result < E{a'^) >— a'^ . 

Now we want to study the Binder's cumulant of a random variable L, where the 
general cumulant ^^271 (in the < L^"+^ >= case which is of interest to us) is defined as 

<L'''^>c 

< >l <L'> 

^4 = ^^ = ^^-3, (A5) 

< > -15 < L2 >< > +30 < >3 

9^ ^JT^. (A6) 

The "natural" estimator for over a set of N measurements is given by 

54=(^i:i?)(^i:i?)"-3, (A7) 

but as we will see this is not a self-averaging quantity, in the sense that < g^ >^ g4, and 
we have to redefine it. To calculate the bias of g^ we introduce the 0{^) fluctuations 772 
and 774 defined by 

i 

^4 = ^E4'-<^'> (A8) 

i 



and we expand Eq. ( |A7| ) around rj = 0, to obtain 



5'4 = 94 



<L^> <L^> < L2 >2 < L2 >< > 



(A9) 



Now to evaluate the corrections we take the expectation value, and we get 



< 5-4 >= 94 



< L2 >2 



< V2 > 



< L2 >< > 



< V2V4 > +^ ( ^ 



(AlO) 



where 



< V2V4 > = - {< > - < X L' >) 



.5=0 



N 



C{L\L\6) + C{L\L\6) 



v<5=0 



2C(L^L^5) - 1 



(All) 



and 



<(x-<x>)(y,+,-<y>)> 
<(x-<x>)(r-<F>)> ' 



(A12) 



If N is much bigger than the correlation times of interest then C(X, F, 5) ~ for 5 > 5, 
with 5 < 77,, so we can write approximately 



N 

<vl> ^1^[<L'>-<L' >') {2tUL\ L') - 1) 

where we have defined the integrated autocorrelation time 

00 

n^t{X,Y) = Y.CiX,Y,5). 

5=0 



< 7]2V4 > ^-{<L'>-<L^><L'>) {tUL\ L") + tUL\ L')-l) (A13) 



(AM) 



(A15) 



In order to correct for the 0{j^) bias we note that in Eq. (|A9|) the terms linear in the 
fluctuations do not contribute to the expectation value. So it is sufficient to define an 
"improved" estimator g\ as 



9 a = 9a 



1-3 



+ 2 



< L2 >2 ' < L2 >< > 



9i 



1-2^^+ 



< L2 > <L^> 



. (A16) 



To the order we are working we can substitute the quadratic fluctuations terms with the 
"naive" estimator for their expectation value, and we obtain 



92 



1 - 



J_y- t2t2 " I J_V /'27'2 

2 



1 E 



2rint(i:',i:')-i 



+ 



+ 



N 



T,^,{L\L') + n^,{L\L'')-l 



(A17) 



To evaluate the variance of we can write using Eq. ( |A16D 



--< (giy > -9l 



9l 



<m> , <vt> 



We can estimate this quantity to the lowest order with 



< L2 >< > 



(A18) 



2 r 



E 




1 E 



+ 



+ 



1 t2tA 



2ri„t(L^L^)-l 



(A19) 



APPENDIX B: DENSITY OF STATES METHOD 

To extract the value of an observable O from a monte carlo simulation we average 
over the sequence of values Oj generated by the algorithm. This is nothing else than an 
approximation of the path integral formula 

Using the fact that we know the probability distribution for the fundamental fields we 
can extrapolate the observable's value around the f] value we have used in the simulation. 
To this effect we note that 

JVijO [ij] exp {-A/3S [tP]) exp {-(33 [V^]) < Oexp {-ApS) > 



<0>il3 + Al3) 
<0>i(3 + Ap) 



JVtPexp {-Af3S [ip]) exp {-(3S [t/;]) < exp {-A/3S) > 

E^0,exp{-A/3S,) 



E^exp{-A/3Si) ' 

so it is sufficient to measure the action Si correlated with each Oi value, and resum each 
measure according to Eq. ( p3|) . The "resummed" value will be reliable only for small 
A/3, because we have a good statistical sampling of action distribution only in a small 
range around the mean value, so we must do a preliminary determination of the allowed 
resummation range. During our simulation we have used essentially the same method as 



Alves et al. in [^]: we plot histograms of the energy distribution, we set a minimum 
and maximum energy by imposing that at least 2.5% of the energy distribution is present 
below the minimum and above the maximum, and we translate this energy range in a /3 
range by finding the (3 shifts that via resummation give the found energy limits. 

If the resumming ranges of different simulations overlap we can "patch" the values 
coming from different extrapolations to obtain a more accurate estimate. 

The determination of an unbiased estimator for the resummed Binder cumulant and 
its variance can be obtained with the same method as in Appendix 0. We only note that 
there are some complications when we combine together different extrapolations, because 
in this case to cancel the bias we must use an improved estimator for the variances, 
depending on many-point correlations which can be estimated reliably only with very 
high statistics. 
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TABLE III. Run parameters (continued) 



AT 


AT 


P 


]\T 

'meas 


Tint 


Rew. range 


a 
O 


5 


l.YoUU 


1do84 


0.7 


1.727 - 


1.772 


o 
O 


p- 

5 


l.YbUU 


40960 


T A 

7.4 


1.738 - 


1.785 


o 



p- 
5 


l.TooU 


40960 


O A 

8.4 


1.741 - 


1.791 


Q 
O 


r. 



1 77nn 


/inQRn 

4UyDU 


Q n 
y.u 


1.745 - 


1.796 


O 




rr 



1. / ( oU 


40960 


inn 
10.0 


1.750 - 


1.801 


8 


5 


1.7800 


57344 


11.3 


1.755 - 


1.806 


8 


5 


1.7850 


40960 


12.0 


1.759 - 


1.811 


8 


5 


1.7900 


40960 


12.4 


1.762 - 


1.815 


8 


5 


1.7950 


40960 


13.1 


1.770 - 


1.822 


12 


5 


1.7600 


39936 


14.0 


1.746 - 


1.774 


12 


5 


1.7650 


39936 


18.2 


1.751 - 


1.777 


12 


5 


1.7700 


39936 


22.9 


1.756 - 


1.783 


12 


5 


1.7750 


39936 


28.5 


1.761 - 


1.789 


12 


5 


1.7800 


23552 


27.5 


1.766 - 


1.792 


12 


5 


1.7850 


23552 


43.1 


1.771 - 


1.798 


12 


5 


1.7900 


23552 


55.8 


1.776 - 


1.804 


12 


5 


1.7950 


23552 


51.6 


1.780 - 


1.808 


16 


5 


1.7600 


28672 


24.0 


1.751 - 


1.767 


16 


5 


1.7650 


28672 


24.2 


1.756 - 


1.774 


16 


5 


1.7700 


28672 


35.6 


1.761 - 


1.779 


16 


5 


1.7750 


28672 


46.4 


1.768 - 


1.784 


20 


5 


1.7670 


18432 


52.4 


1.761 - 


1.773 


20 


5 


1.7690 


18432 


58.5 


1.763 - 


1.774 



TABLE IV. Run parameters (continued) 









-^meas 


Tint 


Rew. range 


8 


6 


1.7750 


32768 


3.9 


1.754 - 


1.799 


8 


6 


1.8000 


32768 


5.4 


1.7775 - 


1.824 


8 


6 


1.8250 


32768 


7.5 


1 son 


1 .ooi 


8 


6 


1.8750 


32768 


11.2 


1.850 - 


1.902 


12 


6 


1.7750 


49152 


5.6 


1.762 - 


1.7875 


12 


6 


1.8000 


49152 


10.1 


1.787 - 


1.812 


12 


6 


1.8125 


16384 


14.7 


1.799 - 


1.826 


12 


6 


1.8250 


49152 


19.5 


1.811 - 


1.837 


12 


6 


1.8375 


16384 


30.6 


1.824 - 


1.850 


12 


6 


1.8750 


49152 


87.9 


1.861 - 


1.889 


16 


6 


1.7750 


28672 


7.5 


1.767 - 


1.783 


16 


6 


1.8000 


28672 


12.0 


1.792 - 


1.809 


16 


6 


1.8125 


16384 


21.9 


1.805 - 


1.820 


16 


6 


1.8250 


28672 


30.1 


1.817 - 


1.832 


16 


6 


1.8375 


16384 


51.2 


1.829 - 


1.846 


16 


6 


1.8750 


28672 


466.3 


1.886 - 


1.884 


20 


6 


1.8200 


16384 


30.3 


1.814 - 


1.825 


20 


6 


1.8250 


18920 


55.0 


1.819 - 


1.831 


20 


6 


1.8300 


16384 


57.6 


1.825 - 


1.836 


24 


6 


1.8220 


12288 


28.5 


1.818 - 


1.826 


24 


6 


1.8250 


27008 


62.9 


1.821 - 


1.829 


24 


6 


1.8280 


27008 


80.2 


1.823 - 


1.832 



TABLE V. Run parameters (continued) 











3 




N 


' lilt 


J. VV • J. C«J.±C^^^ 






7 




1 87nn 




Q91 P> 


O L.O 


l.OUl - l.OOU 






7 








1 1 fifin 


97 9 


1 R(\7 1 RRA 

±.OU 1 - X.OOt: 






7 




1 ssnn 

± .OOVJU 




1 99J^8 


9^^ ^ 


1 ^79 1 


14 




7 




1.8850 




9216 


32.0 


1 878 - 1 892 






7 




1 879 






4n 7 


1 Sfi7 1 877 






7 
t 




1 .0 / (JU 




1 791 


47 R 


1 R7n ^ 

l.O 1 U - l.OOU 






7 




1 S77^ 

1 .0 / 1 o 




1 84nn 


oo.o 


± .O 1 X - X .oo^ 






7 
t 




-L .OOUU 




1 4790 


77 4 


1 R7zL 1 

A..O - l.OOO 






Q 
o 




-L .yUUU 




1 ^^1 9 


1 Q n 


1 SQl 1 QOQ 

± .Oc7l - ± .yuy 


16 




8 




1.9150 




11536 


34.0 


1.907 - 1.922 


16 




8 




1.9200 




13312 


46.2 


1.911 - 1.929 


16 




8 




1.9250 




16384 


39.2 


1.918 - 1.932 


24 




8 




1.9000 




18432 


33.1 


1.896 - 1.904 


24 




8 




1.9150 




29184 


54.5 


1.911 - 1.919 


24 




8 




1.9180 




10144 


27.4 


1.914 - 1.922 


24 




8 




1.9200 




20736 


54.9 


1.916 - 1.924 


24 




8 




1.9220 




10144 


41.4 


1.918 - 1.926 


24 




8 




1.9250 




6144 


54.9 


1.921 - 1.929 


TABLE VI. Run parameters (continued) 


Na 


Nr 






-^mcas 


Tint 








12 


12 




1.375000 


384 


1.0 




0.51106(11) 


0.25133(14) 


12 


12 




1.385000 


384 


1.0 




0.51488(13) 


0.25561(20) 


14 


14 




1.590000 


384 


1.0 




0.596484(95) 


0.35759(14) 


14 


14 




1.595000 


384 


1.0 




0.598358(84) 


0.36011(12) 


14 


14 




1.600000 


384 


1.0 




0.600174(94) 


0.36261(15) 


16 


16 




1.690000 


384 


1.0 




0.631531(63) 


0.40691(11) 


16 


16 




1.695000 


384 


1.1 




0.633037(61) 


0.409088(95) 


16 


16 




1.700000 


384 


1.0 




0.634699(72) 


0.41150(12) 


16 


16 




1.705000 


384 


1.0 




0.636203(78) 


0.41371(10) 


20 


20 




1.769000 


2304 


1.0 




0.654686(17) 


0.440752(27) 


20 


20 




1.770000 


2304 


1.0 




0.654941(13) 


0.441124(22) 


20 


20 




1.827000 


2048 


1.0 




0.669217(16) 


0.462206(26) 


20 


20 




1.828500 


3072 


1.0 




0.669557(10) 


0.462711(17) 


20 


20 




1.830000 


2048 


1.0 




0.669923(13) 


0.463250(22) 


20 


20 




1.872500 


1536 


1.2 




0.679420(18) 


0.477325(29) 


20 


20 




1.877500 


1536 


1.0 




0.680491(14) 


0.478897(24) 


20 


20 




1.915000 


1536 


1.0 




0.688180(16) 


0.490325(27) 


20 


20 




1.925000 


1536 


1.0 




0.690137(15) 


0.493229(24) 



TABLE VIL Run parameters on symmetric lattices, and expectation values of plaquettes 





Na 


N^' 




2 


4 


8 


1.380(4) 


2 


4 


12 


1.380(2) 


2 


8 


12 


1.380(4) 


3 


6 


10 


1.595(1) 


3 


6 


14 


1.5955(8) 


3 


6 


18 


1.5958(5) 


3 


10 


14 


1.5956(12) 


3 


10 


18 


1.5960(7) 


3 


li 


18 


1.5903(15) 


4 


8 


12 


1.6984(7) 


4 


8 


16 


1.6987(5) 


4 


12 


16 


1.699(1) 


5 


8 


12 


1.7669(12) 


5 


8 


16 


1.7678(8) 


5 


8 


20 


1.7682(7) 


5 


12 


16 


1.7684(15) 


5 


12 


20 


1.769(12) 


5 


16 


20 


1.769(3) 


6 


8 


12 


1.812(4) 


6 


8 


16 


1.817(2) 


6 


8 


20 


1.822(2) 


6 


8 


24 


1.824(1) 


6 


12 


16 


1.823(5) 


6 


12 


20 


1.828(2) 


6 


12 


24 


1.825(1) 


6 


16 


20 


1.830(4) 


6 


16 


24 


1.825(2) 


6 


20 


24 


1.824(3) 


7 


14 


20 


1.8747(20) 


8 


16 


24 


1.920(4) 


TABLE VIII. Crossings for 54 curves 








Te/Ai 


2 




1.380(4) 


8.81(8) 


3 




1.59624(13) 


9.889(3) 


4 




1.699(1) 


9.526(30) 


5 




1.76948(3) 


9.0562(7) 


6 




1.8287(11) 


8.73(3) 


7 




1.8747(30) 


8.38(6) 


8 




1.920(5) 


8.2(1) 



TABLE IX. Critical temperature by two loop asymptotic scaling formula 



Nr 






A ri 




3 


1.588(2) 


9.70(5) 


0.76(4) 0.35(2) 


5.4 


4 


1.687(2) 


9.25(3) 


0.66(4) 0.40(2) 


15.3 


5 


1.759(3) 


8.82(5) 


0.47(2) 0.38(2) 


4.0 


TABLE X. Results of the fit 




Pi 




Pl 


Ai/Al 


2 


1.380(4) 




1.8800(30) 


3.599(46) 


3 


1.59624(13) 


2.1768(30) 


4.470(35) 


4 


1.699(1) 




2.2986(6) 


4.713(14) 


5 


1.76948(3) 




2.3726(45) 


4.766(56) 


6 


1.8287(11) 




2.4265(30) 


4.709(39) 


8 


1.920(5) 




2.5115(40) 


4.644(77) 


TABLE XL A ratios from asymptotic formula 


Nr 






T/Al 


T/A^ 


2 


0.1693(3) 


0.2466(6) 


29.7(2) 


2.27(2) 


3 


0.1462(2) 


0.2039(3) 


41.4(3) 


2.89(2) 


4 


0.13848(4) 


0.19022(7; 


) 42.1(1) 


2.848(4) 


5 


0.1342(2) 


0.1828(3) 


40.6(5) 


2.69(2) 


6 


0.1312(2) 


0.1776(3) 


38.7(3) 


2.54(2) 


8 


0.1267(2) 


0.1700(3) 


36.0(4) 


2.32(2) 


16 


0.1162(4) 


0.1526(6) 


32.0(8) 


1.97(4) 




TABLE XII. Perturbative 


,: Wilson case 




Nr 






T/Ai 


T/A^ 


2 


0.2306(7) 


0.286(1) 


8.81(8) 


1.49(1) 


3 


0.19941(2) 


0.24085(3) 


9.889(3) 


1.627(6) 


4 


0.1873(1) 


0.2239(1) 


9.526(30) 


1.549(2) 


5 


0.179889(3) 


0.213616(4) 


9.0562(7) 


1.461(1) 


6 


0.1741(1) 


0.2057(1) 


8.73(3) 


1.398(2) 


7 


0.1698(3) 


0.1998(4) 


8.38(6) 


1.34(1) 


8 


0.1658(4) 


0.1944(5) 


8.2(1) 


1.30(1) 




TABLE XIII. Perturbative «„: 


Symanzik case 









Ti/Ai 


/?2 


T2/A2 


1.380(4) 


1.1280(35) 


1.832(15) 


1.331(30) 


2.953(21) 


1.59624(13) 


1.36942(16) 


2.15542(83) 


1.55480(16) 


3.3651(13) 


1.699(1) 


1.5026(13) 


2.2247(73) 


1.6876(13) 


3.485(12) 


1.76948(3) 


1.59156(3) 


2.2072(2) 


1.77784(3) 


3.4781(3) 


1.8287(11) 


1.6628(13) 


2.1876(70) 


1.850(1) 


3.462(11) 


1.8747(30) 


1.7163(34) 


2.136(18) 


1.9041(34) 


3.391(29) 


1.920(5) 


1.7674(56) 


2.119(29) 


1.9557(56) 


3.371(47) 



TABLE XIV. Data for the effective schemes 




FIG. 1. Binder cumulant: lattices with Nr = 2 




FIG. 2. Binder cumulant: 



lattices with Nr = 3 




FIG. 3. Binder cumulant: lattices with Nr = 4, 




FIG. 5. Binder cumulant: lattices with AT^ = 6 
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FIG. 6. Binder cumulant: lattices with Nr = 7 




FIG. 7. Binder cumulant: lattices with Nr = 8 



X 



o 3 X [6, 10, 14, 18] 
□ 4x [8, 12, 16] 
A 5 X [8, 12, 16, 20] 
• 6 X [12, 16, 20, 24] 
■ 7 X [14, 20] 
A 8 X [16, 24] 
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FIG. 8. Universal scaling behavior of Binder cumulant 
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FIG. 9. Universal scaling behavior of Polyakov loop 
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FIG. 10. Plot of the Tc/A ratio by use of the "bare" scheme 




FIG. 11. Critical behavior of Polyakov line 




FIG. 12. "Bare" scheme corrected with the first sub-asymptotic coefficient of (3 function 




FIG. 13. Comparison of "bare" and effective schemes 



